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We investigate the physical meaning of quantum superposition states between reactants and 
products in electron transfer reactions. We show that such superpositions are strongly suppressed 
and to leading orders of perturbation theory do not pertain in electron transfer reactions. This 
is because of the intermediate manifold of states separating the reactants from the products. We 
provide an intuitive description of these considerations with Feynman diagrams. We also discuss 
the relation of such quantum coherences to understanding the fundamental quantum dynamics of 
spin-selective radical-ion-pair reactions. 



I. INTRODUCTION 

The role of quantum superpositions in irreversible 
chemical reactions, in particular electron transfer reac- 
tions [H, has received little attention. The specific ques- 
tion we will here address is this: to first-order reaction 
kinetics, can there be quantum superpositions between 
reactants and reaction products of an electron transfer 
reaction? Most if not all electron transfer reactions are 
known to occur in two steps, as shown in FigfTJ depicting 
the potential curves of the reactants and the products: 

(1) the electron tunnels from the initial state \R) into an 
intermediate state \P*), which is a vibrational excited 
state of the product, the ground state of which is \P). 

(2) The excited state \P*) decays to the ground state 
\P), with the emission of a photon or a phonon. This 
decay usually takes place at a rate much faster than the 
first step, so it is the latter that mostly determines the 
total reaction rate. Most works on electron transfer deal 
with calculating the rate of this first step [1-0] • More- 
over, relevant considerations [8] of quantum coherence 
have considered, again, this first step, and of course have 
appeared in the much broader context of chemical reac- 
tions [9\. To our knowledge, however, the specific ques- 
tion we defined, i.e. whether the reactant state can be 
coherently coupled to the final product state in electron 
transfer reactions has not been addressed so far. 
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FIG. 1: (a) Level structure of reactants and reaction products 
in the electron transfer taking place in charge recombination 
of radical-ion pairs (reactants) and the formation of neutral 
products. 



This question can be rephrased: Can states of the form 
\ip) = Ci|R) + C2IP) be physically realized to first-order 
reaction kinetics? The intuitive answer is negative, be- 
cause of the existence of the intermediate states IP*): a 
substantial amplitude for the product state \P) would 
translate to an equally suppressed amplitude of the in- 
termediate \P*) and thus a negligible amplitude of the 
reactant state \R) . Here we are going to prove that the in- 
tuitive answer is indeed correct. In retrospect, this might 
seem to be an obvious statement. 

However, it can't be that obvious after all, for the fol- 
lowing reasons. The seemingly innocuous question we 
will here address is directly related to understanding 
the fundamental quantum dynamics of spin-dependent 
recombination reactions of radical- ion pairs [lOj . which 
have recently received a lot of attention from the quan- 
tum information science perspective [lll - [l9j . In par- 
ticular, the traditional theoretical description of spin- 
selective radical-ion-pair (RP) reactions has been ques- 
tioned, and new approaches have been put forward by 
the authors in (Til \l2L llil ] , who introduced quantum mea- 
surements as conceptually important for understanding 
the quantum dynamics of these reactions. These papers 
have led to an ongoing debate [HI IT7I I20I - I23I ] which at the 
least illuminates the richness of the topic in terms of the 
quantum mechanical concepts involved. Currently there 
exist at least three different theories, all purported to be 
derived from first principles, attempting to describe RP 
reactions. The traditional theory, sometimes called the 
Haberkorn master equation [24j], was explicitly derived, 
for the first time to our knowledge, by Ivanov et al. [22j , 
and supported by Purtov along similar lines of thought 
23]. Both of these papers, however, start out from this 
seemingly innocent assumption: the existence of quan- 
tum superpositions between the reactants (the radical- 
ion-pair state) and the reaction products (the neutral 
product state) of an electron transfer reaction. Further- 
more, the considerations in also start out from the 
assumption that the radical-ion-pair state can be coher- 
ently coupled to the neutral product state. Thus, it is 
fair to say that what we are going to prove here is not an 
obvious, but a rather subtle physical statement regard- 
ing electron-transfer reactions in general, and radical-ion- 
pair reactions in particular. 

In Section II we will rigorously prove our main as- 



sertion, namely that to first order reaction kinetics, the 
product state of an electron transfer reaction cannot be 
coherently coupled to the reactant state, exactly due to 
the existence of intermediate states. In Section III we 
will elaborate on the physical consequences of neglecting 
these states, or stated equivalently, treating RP reactions 
by the Haberkorn master equation. We will show that 
according to the latter, the evolution of the entropy of 
the RP spin state exhibits an unphysical behaviour, not 
evidenced by the other two theories based on quantum 
measurement considerations. 



II. REACTANT-PRODUCT QUANTUM 
COHERENCE IN ELECTRON TRANSFER 
REACTIONS 

We will consider a simplified version of the problem 
keeping all the ingredients that are essential for the ba- 
sic physics. Reducing the problem to an electronic one, 
i.e. neglecting the nuclear degrees of freedom, we depict 
in Fig. [2] the energy levels of the system comprised of 
the single reactant state \R), a manifold of intermediate 
states |P*) to which the electron can tunnel from \R), 
and the single product state |P) to which |P*) can de- 
cay by the emission of a photon or a phonon. Letting 
\P) define the zero of the energy scale, the unperturbed 
Hamiltonian of the system is 



R Q = ujRa' R a R 



Ljp*a P ,ap* 



^Wkaj^ak (1) 



The sum over k denotes all modes of the photons 
(phonons) with frequency Wk = uir due to energy con- 
servation. We will use fermionic operators to denote the 
creation (a^) or annihilation (a) of a single occupation of 
the relevant states. The eigenstates of R are denoted 
by \R, P*,P, iVk), where R, P*,P — 0, 1 are occupations 
of the reactant, the intermediate and the product states, 
respectively, and similarly iVk = 0, 1 denotes the occupa- 
tion of the photon state with wave vector k. For example, 
the initial state is \i) = |1; 0; 0; 0), whereas the final state 
is one among the states |0;0;l;lk), where k runs over 
all photon wave vectors satisfying energy conservation. 
These states will be collectively denoted by |P), i.e. with- 
out reference to the index k. The tunneling Hamiltonian 
coupling \R) to |P*) is simply Rt — J2i ^i a R a p* + c - c - 
Finally, the decay Hamiltonian coupling |P*) to the ra- 
diation (or phonon) reservoir and leading to the ground 
state \P) is Hd = J2i J2k Ci^ap'apa^ + c.c. The coupling 
constants c^k readily result from the quantum theory of 
light-matter interaction, and are c^k = \/hid^/2VA4i, 
where V is the volume in which the radiation field modes 
are defined, and Aii the relevant matrix element for 
the transition P* — > P. The total interaction Hamil- 
tonian is Vi n t = Rt + Rd- It is worth considering two 
simple cases before proceeding: (i) Neglecting the pres- 
ence of |P) and considering only the tunneling transition 
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FIG. 2: Simplified level structure of an electron transfer pro- 
cess. The tunneling Hamiltonian Rt couples the "reactants" 
state \R) with the manifold of excited states |P*), which then 
radiatively decay to the ground state |P), the "product" state. 



\R) —> \P*), moreover assuming that A.; = A for all i, and 
using Fermi's golden rule, it is found within l st -order 
perturbation theory that this transition's rate is k = 
(2Tr/h 2 )\\\ 2 p*(ujp), where p*{ui) is the density of states of 
the {P*}-manifold calculated at the energy of the reac- 
tants ujr. (ii) Neglecting the presence of |P) and consid- 
ering only the radiative transition |P*) — » |P), moreover 
assuming that Mi = M for all i, and again using Fermi's 
golden rule, it is found within l st -order perturbation the- 
ory that this transition's rate is T = (2ir / 'fi, 2 )|.M| 2 p(ui p*) , 
where p(ui) is the density of states of the radiation field 
calculated at the energy of the excited intermediates up* . 
It is noted that in realistic systems it usually is V ^S> k. 

We will now move to consider the action of Vi n t on the 
entire system using time-dependent perturbation theory. 
The initial state \R) = |1; 0; 0; 0) will evolve to |V>t) = 
c R (t)\R) + I>p,(i)l^*) + c P (t)\P). To th -order per- 



turbation theory it is c"p (t) = 1 and c 



(0), 
F" 



(t)=cp(t) = 0. 



It is clear that within l st -order perturbation theory, Vint 
can couple |P) to |P*) through Rt, but not to |P), hence 
cp 1 "* = 0. Let's now move to 2 nd -order perturbation the- 
ory, where a possible coupling between |P) and |P) ap- 
pears through the product RdRr- Assuming that the 
reactant and intermediate states are non-degenerate, i.e. 
uip ^= up* it follows from 2 nd -order time-dependent per- 
turbation theory (setting 5^ = lu^ — up) that 
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As the manifold of intermediate states is a dense quasi- 
continuum, and as mentioned before, Xi — A and c^.k = 
Ck for all i can be taken out of the sum, moreover since 
l/(wp* — uir) is an odd function of cjp» — ujr, it follows 

(2) . 1 

that c P — 0. Even in the case of resonant tunneling, i.e. 
supposing the existence of one intermediate state |P*) 
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resonant with \R), it follows that 

» iSkt (l - iS k t) - 1 
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where the second equality follows in the limit <5k — > 
imposed by energy conservation. The phase e l<5kt will 
average c p to zero. Even if 5k = identically, this 
contribution is not first-order in the reaction kinetics, 
but third, so it can be neglected. 

Now it is also clear that coherently coupling | P) to | R) 
is possible only within even-order perturbation theory, 
hence the first non-zero term in the sought after quantum 
superposition appears at 4 th -order perturbation theory, 
easily found to be negligible. There is a very transpar- 
ent way to summarize the above results with Feynman 
diagrams, shown in FigJ3J Not taking into account the 
intermediate states |P*), i.e. supposing a physically ab- 
sent direct coherent coupling between \R) and \P), is 
depicted by the diagram shown in Fig|3^,. Taking the 
full level dynamics into account, we have shown that the 
2 nd -order perturbation theory diagram depicted in Fig(3jD 
has a zero amplitude due to the energy-odd denomina- 
tor. The first non-zero amplitude for coherently coupling 
\R) to \P) results from the 4 th -order diagram of Fig(3J;. 
In this diagram, the intermediate state |P*) emits a vir- 
tual photon that is re-absorbed by the produced product 
state, leading again to the intermediate manifold before 
the actual product is generated. Finally, the electron 
transfer reaction actually proceeds within first order re- 
action kinetics through two consecutive real transitions, 
i.e. there is a reaHransition from \R) to IP*), followed by 
another real transition from \P*) to \P). We symbolize 
this with the diagram shown in Fig[3Ji, where the dashed 
line crossing the \P*} states is supposed to mean that the 
\P*) states are really, as opposed to virtually, populated 
by the real transitions from \R). 



III. NEGLECTING INTERMEDIATE STATES 

As mentioned in the introduction, the authors in [22[ 
and [23l ] derived the traditional master equation of spin 
chemistry treating the recombination process according 
to the diagram in Fig. [3^, i.e. by neglecting the interme- 
diate states \P*). This is an example of what field theory 
calls an " effective theory" . There are many examples in 
physics, especially in nuclear physics [25j . where effec- 
tive theories are indeed a good approximation, destined 
however, to fail in some regime of the relevant parameter 
space. For example, weak interactions were initially de- 
scribed by the effective Fermi theory, according to which 
the neutron decays as shown in FigH^,. After the de- 
velopment of electroweak theory, the relevant Feynman 
diagram was found to include the intermediate W~ bo- 
son (Fig|l]b). The Fermi theory was a good low energy 
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FIG. 3: (a) "Effective" description of electron-transfer reac- 
tions neglecting the intermediate states and assuming a direct 
coupling between reactants and products, (b) Taking into ac- 
count the intermediate states, this is the lowest-order (2 nd ) 
coherent coupling between reactants and products, the ampli- 
tude of which is zero, (c) The first non-zero amplitude comes 
from this 4 th -order diagram, (d) Within frist-order reaction 
kinetics, the reaction proceeds by a real transition from \R) 
to \P*) (indicated by the dashed line), followed by another 
real transition from \P*} to \P). 
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FIG. 4: (a) Early understanding of neutron decay, leading to 
several problems, such as an infinite electron-neutrino scat- 
tering cross section, (b) The infinity was remedied when the 
actual vertex involving the intermediate W~ boson was un- 
derstood from electroweak theory. 



approximation, plagued, however, by several problems, 
such as an electron- neutrino scattering cross section tend- 
ing to infinity at high energies. This and other problems 
were remedied by the full electroweak theory. Similarly, 
the traditional master equation of spin chemistry is a 
very good approximation, plagued however, with several 
conceptual problems. We will now turn to those. 
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FIG. 5: Comparison of the three master equations by Komi- 
nis (solid line), Jones-Hore (dashed line) and the traditional 
(dotted line) in predicting the von-Neuman entropy S(p). The 
density matrix p was normalized by Tr{p} in order to remove 
the artefact of diminishing RP population due to recombina- 
tion. 



locked in the non-reactive triplet state forever (we took 
kr = 0). Indeed this is what the prediction of all three 
theories is, albeit with a quantitative difference of what 
exactly this fraction is [l2j . This difference is immaterial 
for the current consideration. What is also expected on 
physical grounds is the fact that while at the beginning 
and the end the state of the radical-ion pairs is pure, i.e. 
at t — we have a pure singlet state and at t — >• oo we 
have a pure triplet, during intermediate times we have 
a mixture of pure states. Hence one would expect that 
the von-Neuman entropy S(p) = — Trjplnp}, where p is 
the density matrix describing the two-electron spin state, 
would start from zero, acquire a non-zero value at inter- 
mediate times, and return to zero. In FigjSJi we compare 
the predictions of the three master equations. It is clear 
that the two theories based on the quantum measure- 
ment approach to the recombination of RPs predict the 
above described behavior expected on general physical 
grounds. In contrast, S(p) = at all times according to 
the traditional theory. 



We consider a very simple example of a radical-pair 
consisting of just two unpaired electrons and no nuclear 
spins. We furthermore set to zero all magnetic interac- 
tions and consider a single recombination channel, e.g. 
the singlet, with ks being the recombination rate. We 
take as an initial state the maximally coherent state 
\ip ) = (\S) + \T ))/V2, where IS") and \T ) are the two- 
electron singlet and zero-projection triplet state, respec- 
tively. What is physically obvious is that a fraction of the 
radical-ion pairs, initially in the state \ipo) will recombine 
through the singlet channel, while the rest will remain 



The aforementioned problem of the Haberkorn master 
equation can be traced to the absence of intermediate 
states in the theoretical paradigm underlying the tradi- 
tional theory. The presence of such states leads to a fun- 
damental and unavoidable S-T dephasing within first or- 
der reaction kinetics through Feynman diagrams like the 
ones presented in [T2| . Concluding, we have explored the 
physical existence of quantum superpositions between re- 
actants and products in electron transfer reactions and 
connected the result to the fundamental quantum dy- 
namics of spin chemistry reactions. 
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